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Abstract 



We present a new convergent iterative solution for the two lowest quantum wave 
functions ip ev and t/w of the Hamiltonian with a quartic double well potential V in 
one dimension. By starting from a trial function, which is by itself the exact lowest 
even or odd eigenstate of a different Hamiltonian with a modified potential V + SV, 
we construct the Green's function for the modified potential. The true wave func- 
tions, tp ev or ip d, then satisfies a linear inhomogeneous integral equation, in which 
the inhomogeneous term is the trial function, and the kernel is the product of the 
Green's function times the sum of SV, the potential difference, and the correspond- 
ing energy shift. By iterating this equation we obtain successive approximations to 
the true wave function; furthermore, the approximate energy shift is also adjusted 
at each iteration so that the approximate wave function is well behaved everywhere. 
We are able to prove that this iterative procedure converges for both the energy 
and the wave function at all x. 

The effectiveness of this iterative process clearly depends on how good the trial 
function is, or equivalently, how small the potential difference SV is. Although each 
iteration brings a correction smaller than the previous one by a factor proportional 
to the parameter that characterizes the smallness of SV, it is not a power series 
expansion in the parameter. The exact tunneling information of the modified po- 
tential is, of course, contained in the Green's function; by adjusting the kernel of 
the integral equation via the energy shift at each iteration, we bring enough of this 
information into the calculation so that each approximate wave function is expo- 
nentially tuned. This is the underlying reason why the present method converges, 
while the usual power series expansion does not. 



PACS: 03.65. Ge, 11.15. Me 
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1. Introduction 



The quartic potential in one dimension with degenerate minima, 



V(x) = ^ V - a 2 ) 2 



(1.1) 



is the archetypical bound-state tunneling problem in quantum mechanics. Its perturbation se- 
ries is well known but does not converge as one would expect since the exponential contribution 
from tunneling between the two minima is necessarily nonperturbative. 

An alternative form of the same problem is obtained by setting q = y/2ga(a — x) so that 
the Hamiltonian becomes 



This shows that the dimensionless (small) expansion parameter is related to l/y/8 ga 3 ; as 
it turns out, the relevant parameter is its square. In this paper we shall take a = 1 so that the 
expansion parameter is 1/g; in the literature one often finds the assumption 2ga = 1 (placing 
the second minimum of the potential at q — 1/g) so that 1 / ^8ga 3 reduces to g and the 
anharmonic potential appears as (l/2)g 2 (l — gq) 2 . Then g appears with positive exponents 
instead of negative, but the coefficients of the power series are the same as with our form of 
the potential, apart from the overall factor 2ga. 

Recent work[l-10] on this problem has concentrated on the Feynman path integral by which 
the nonperturbative tunneling effect is represented by a series of instantons and anti-instantons. 
Emphasis is sometimes on understanding the asymptotic behavior of the perturbation series as 
a signature of the tunneling influence, and sometimes on supplying a nonperturbative correction 
to supplement the truncated perturbation series. 

In particular, Zinn-Justin[4-7] has produced an elegant implicit equation that determines 
all the eigenenergies in terms of two functions given as power series: one series is the Brillouin- 
Wigner expansion for the energy, and the other is the steepest-descent expansion for the action 
of a single instanton. To any order in the power series, an expression for the energy is obtained 



H = -\ + V(x) = 2gah, 



(1.2) 



where 




(1.3) 
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which, through the implicit equation, contains information pertaining to arbitrary numbers 
of instantons. One is thus in the odd position of possessing approximations to finite order in 
/ = l/(ga 3 ) but in some sense correct to all orders in e^ 1 ^ . 

We offer here a completely different method which leads to a result of similar character 
although its equivalence with that of Zinn-Justin is yet to be studied. Our method is based 
entirely on the Schroedinger equation, with no reference to path integrals or to analytic con- 
tinuation in complex g. We start from an approximation to the desired wavefunction, which 
itself satisfies exactly a Schroedinger equation with a modified potential V + SV. From this 
trial wavefunction we obtain the Green's function for the modified potential. The true wave- 
function is then given implicitly by acting with the Green's function on a source obtained 
by multiplying the difference between the two potentials by the true wavefunction itself. By 
iterating this equation we obtain successive approximations to the true wavefunction. We are 
able to prove that these approximations actually converge. 

The effectiveness of this iterative procedure clearly depends on how good the trial function 
is, or equivalently, how small the potential difference 5V between the modified potential and 
the true potential is. Although each iteration brings a correction smaller than the previous 
one by a factor proportional to the parameter that characterizes the smallness of SV, it is not 
a power series expansion of the parameter, since the Green's function also contains the exact 
tunneling information of the modified potential. This is the reason why the present method 
converges, while the usual power series expansion does not. 

Another important feature of the method is that the approximate energy is also adjusted 
at each stage so as to prevent the approximate wavefunction from getting too large at the 
second minimum. Thus we simultaneously develop a convergent sequence of approximations 
for the energy. Like that of Zinn-Justin, the sequence is correct to successively higher orders 
in / (or, 1/g in our notation), but each step is in some sense correct to all orders in e" 1 ^. 

Because we are directly dealing with the Schroedinger equation 

HiP(x) = EiP(x), (1.4) 

the proof of convergence is explicit, and it follows the traditional method of analysis, as we 
shall see. 
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Recall that the formal perturbative expansion of 

tP(x) = e - 9S(x) 



(1.5) 



and its associated energy E may be written as 



and 



gS = gS + S 1 + -S 2 + --- (1.6) 



E = gE + E 1 + -E 2 + ■ ■ ■ , (1.7) 
9 



where (after setting the parameter a = 1 in the potential V(x)) 

So(*) = i(*-l) 2 (z + 2), S 1 (x)=ln^, S 2 {x) = ±- *\\ r ••■(1.8) 
and 

E = l, E l = -- E2 = ~, .... (1.9) 
4 64 

Both expansions (|1.6| ) and ( |1.7| ) are divergent; furthermore, at x — — 1 and for n > 1, 
each 5 n (x) is infinite. The reflection x — > — x gives a corresponding asymptotic expansion 
<Sn(^) — > S n (—x), in which each S' n (— x) is regular at x = —1, but singular at a; = +1. 

In this paper, we will not use the perturbative expansion. Instead, we will follow a new 
approach suggested in Section 3 and 4 of a previous paperfll] ( referred to as I hereafter ). 
Accordingly, we will construct an appropriate trial function to start a new iterative procedure 
to derive the lowest even, or odd, state which satisfies 

Hip ev = E ev ip ev or Hi/j od = E od ip od . 

As first step, we retain the first two terms in the expansion (|1.6| ) of S. For x > 0, the function 

<j) + (x) = e- sSo(l) - Sl(l) 

is well behaved and, for g large, approximates the overall behavior of ip ev . However, + does 
not satisfy the boundary condition of ip ev at x — 0. 
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In Section 2, we use + to construct a trial function 0: 

for x > 



<f>(x) 



4 , 



if < x < 1 

(l + ge-^)0 + (x), if x>l 



where 



x 



) = e -»S (-x)-Si(x) _ e -fs+gSo(:r) 



1 + x 



for rr < 

= (f)(—x). 

Thus, is even in x and at x — 0, its derivative 0'(O) = 0. Furthermore, it can be verified 
that (f>(x) and its first derivative <f)'(x) are continuous at all x, and in addition, (f>(x) is the 
lowest eigenstate of a Schroedinger equation with a different Hamiltonian: 

1 d 2 



where except for x near ±1, 



— — — = O(-) or smaller. 



When x is near ±1, has the correct Gaussian behavior, like ip ev (x). 

Next, we construct the Green's function (x\D\y), whose matrix form is defined by 

D = -2(1)9 (j)- 2 9(f), 

in which the matrix element of 9 is a step function: 

0, for x > y 



(x\9\y) 



-1, for x < y. 



Thus, D satisfies 



1 d 2 



- 2^ + v ( x ) + w ( x ) - a\ ( x \D\y) = K x - y) 



and, in the matrix notation, 



ip ev = <p + D(w + E ev - g)ip e 



The iterative procedure beginning from to arriving at ip ev is discussed in Section 2. The 
convergence of this procedure is proved in Section 3. 

To demonstrate the flexibility of our approach, we give a variation of the procedure for 
the derivation of the lowest odd eigenfunction ip d of the double-well potential problem. In 
Section 4 we begin with the trial function </>+. As we shall see, using <j> + , instead of 0, as the 
zero th approximation, we can arrive at a solution ip+(x) of the differential equation 

H^ + (x)=[-~ + V{x)]Mx) = E + Mx) 
with ip+(x) and its derivative ip' + (x) satisfying the boundary conditions: 

V>+(o) M») 

and 

at x = oo, tp + (oo) = 0. 

(At x = — oo, ip+(x) diverges.) For x > 0, the derivation of tp + from + follows the same 
iterative steps as those in Sections 2 and 3, through the replacement of by + and tp ev by 
Likewise, the iteration is a convergent one. 
In Section 5, we start from ip + to first construct an odd trial function x- The function \ 
is defined to be: 
for x > 

X = $+-il>-, 

where (in the matrix form) 

with the constant 7 determined by the condition: at x — 0, x(0) = 0; 

for x < 

X(x) = -x(~x). 

It will be shown that xi x ) ' s the lowest eigenstate of a new Hamiltonian: 

(H + u)x = E +X , (1.10) 
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where H is the same double-well Hamiltonian, as before. As we shall see, v is 0(e), where 

4 

e = e~3 9 is the very small tunneling parameter, and E + is between the two eigenvalues E ev 
and E od of H. 

Using x as the new zero th order approximation for ip d we first construct the corresponding 
Green's function G, defined by 

g = ~2 X e x - 2 e x . 

Because of (|1.10| ), G and ip Q d satisfy 

(H + u-E + )G= 1 

and 

= X + G{v + E od - E + )^ od . (1.11) 



The integral equation (|1.11| ) now serves as the basis of a new iterative procedure using 



G(u + E od — E + ) as the kernel for iterations. As will be shown, each iteration improves the 
solution towards ip Q d by a factor 0(e). It will be proved that this new iterative sequence is also 
a convergent one. Furthermore, the convergence is point-wise in all x ( as are the previous 
iterations, from (ft to ip ev and from <p + to ip + ). 

Throughout the paper, most of the iterative functions are only known through their func- 
tional relations. To achieve the convergence proof, we rely on deriving the explicit upper and 
lower bounds for a number of these functions and their related parameters. Some of their 
derivations can be rather lengthy; these are given in Appendices A to C. 

From a practical point of view, once a good trial function is constructed, the new iterative 
procedure is quite effective in generating better and better approximations, as will be discussed 
in a particular example given in Appendix D. 
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2. The Iterative Process 



With a = 1, the one-dimensional Hamiltonian ( |1.2| ) becomes 

= T+^V-1) 2 , (2.1) 
where T = — Let be the lowest eigen-state, which is even in x and satisfies 

Hip ev = E ev ip ev . (2.2) 

As mentioned in the introduction, we shall solve ip ev by starting from a trial function cf)(x) 
which satisfies a Schroedinger equation with a different Hamiltonian say, instead of ( |2.2|) , 

(H + w)4> = (E ev + £)(j). (2.3) 

The function 0(x), like ip ev (x), is even in x and goes to at x = ±oo. From 0(x), one can 
readily construct an irregular solution 0(x) which satisfies the same equation (|2~3|): 

0(x) = <p(x) / <fr~ 2 (y)dy 



The standard Green's function (x|D|y) for the Sturm-Liouville type problem[13] can be formed 
by using the product of the bilinear expression 

<p{x)4>(y) - 0(x)0(y) 

times a step-function 9 in x — y. As will be discussed, the eigen-function ip ev of the original 
problem fl2~lp- (|2.2| ) can then be derived as a power series of D. The success of the method 
clearly depends on the choice of the trial function (f)(x). The difference w between these two 
Hamiltonian in (|2.2| ) and ( |2.3|) can be obtained by a straightforward differentiation of 0, as 
we shall see. 

Following the steps outlined in the Introduction, we introduce for x > 0, 

gS (x) = |(x-l) 2 (x + 2), (2.4) 
Si(x) = ln^, (2.5) 

<j>-{x) = e -f s »(-^- s iW = e -|9e + 9 s »W(^-) (2.7) 



and 

f 6Jx) + a=k<J>Jx), for < x < 1 

The functions gSo(x) and Si(x) are identical to those used in I and Ref.[12] (referred to as II 
hereafter; note that the <p + and </>_ in this paper are different from those defined in II.). By 
differentiating + , we find that it satisfies 

(T + V + u)<P + = #+, (2.9) 

where 

V(x) = \g\x 2 - l) 2 (2.10) 

and 

u( X ) = — Ljj. (2.11) 

A good trial function 0(x) for ij) ev (x) must also be even in x. Thus, in the upper expression 
in ( |2~8| ) the additional term proportional to 4>-(x) for < a; < 1 is to make the derivative of 
4>{x) vanish at x = 0; i.e., 

0'(O) = 0. (2.12) 

(Throughout the paper, prime denotes ^, and g is > 1.) At a; = 0, </>+(0) = <f>-(0) = e~i 9 . 
As x increases from 0, d>+(x) increases to 1 at x — 1, while d)-(x) decreases to e~^ 9 at x — 1. 



We modify 0(x) to the lower expression in ( ^8|) so that 0(x) becomes proportional to <p+(x) 
for x > 1; furthermore, 4>{x) and </>'(a;) are both continuous at x — 1. 
The function satisfies 

(T + 1/ + u>)0 = gcf), (2.13) 

where, for x > 0, 

w(x) = u(x) + (2-14) 
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with u(x) and g(x) both even in x. For x > 0, u(x) is given by (|2.11| ) and 

f 2q (g-De 29SoW -^ a for o < x < 1 

£(*) = ^ (s+ i) +M e^<«H»' - (2.15) 

[ for a; > 1, 

which is discontinuous at x = 1. For a; < 0, w(x) = w(—x). 
Introduce the matrix 9 whose matrix elements are 

Following the Green's function notations in I, we define the matrices 

D = -29(f)- 2 e(j) 2 (2.17) 

and 

D = fiD <j)- 1 . (2.18) 
Differentiating D with respect to x from the left, we find 

D' = -2cp'6(j)- 2 6(j) + 2(j)- 1 0(j> 

and 

D" = -2<j)"6(j)- 2 6<p - 2. 

Thus, 

(T + V + w - g)D = I = unit matrix, (2.19) 



and we see that ip ev of (|2T2j ) satisfies 

7p ev = (j) + D(w - £)^ ev , (2.20) 

with its eigenvalue given by 

E ev = g-£. (2.21) 
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In addition, 

oo 

/ w4>ip ev dx 

8 = ^ . (2.22) 



The proof is identical to that given for (4.92) of I. In the following, we give an alternative 
proof by introducing / through 

Ipev = <t>f- (2.23) 

Thus, / satisfies 

f = l+D(w-£)f; (2.24) 
i.e., on account of (|2.16| ) and ( |2.17|) , 



oo oo 

2/ \ r , f A2l 



f(x) = 1-2 J r\y)dy J P(z)(w(z) - S)f(z)dz. (2.25) 

x y 

Since at x = 0, ^(0) = and 0'(O) = 0, we also have /'(0) = 0, which leads from ( gggp to 

(j) 2 (x)(w(x) — £)f(x)dx = 0, 



o 



and gives ( |2.22| ) 



In what follows, we shall solve / and 8 by the iterative sequences {/„} and {£ n }'- 

fo = 1, fl, f2, f3, ■ ■ ■ , fn, ••• 

and 

S = 0, S u £2, £3, £n, ■■■ (2.26) 

We require 

f n = l+V(w -£ n )fn-l (2.27) 

with 

£n = [wf n ^]/[f n ^\, (2.28) 
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where [F] of any function F(x) is defined to be 

oo 

[F] = J (j) 2 (x)F(x)dx. (2.29) 
o 

Because of (P2BD , at x = 

£(0) = 0. (2.30) 

In explicit form, 

oo oo 

h{x) = 1-2 J r 2 (y)dy J (f> 2 (z)(w(z) - S 1 )dz, (2.31) 

x y 
oo oo 

£ x = J <p 2 (x)w(x)dx I J 2 (x)dx, (2.32) 




oo oo 
-2/„,\ j„, /" j2 



/ 2 (a;) = 1-2 J r'(y)dy J r(z)(w(z)-£ 2 )f 1 (z)dz, (2.33) 

x J/ 

oo oo 

£ 2 = J (p 2 (x)f 1 (x)w(x)dx^ J (p 2 (x)f 1 (x)dx, (2.34) 



etc. Because of (|2.28| ), (|2.27| ) can be written either as 

oo oo 

fn(x) = 1 - 2 / <p-\y)dy f 4> 2 (z)(w(z) - S n )f n ^(z)dz, (2.35) 
or in an equivalent form 

fn(x) = / n (0) - 2 y 0- 2 (y)rfy J <P' 2 (z)(w(z) - £ n )f n _ x {z)dz (2.36) 



with 

oo oo 

E n = J 4> 2 (x)f n ^i(x)w(x)dx I J 4> 2 (x)f n ^i(x)dx. (2.37) 
o o 

As we shall show, this iterative procedure is a convergent one. 



x y 
-2/„\j„ f iii 
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3. Some General Properties 



For clarity, we list the following properties of /„ and 8 n in the form of several theorems. 
Theorem 3.1 At all x > 0, the derivative of each f n (x) is negative; i.e., 

fjx) < 0. (3.1) 

Thus, 

/n(0) > f n (x) > f n (00) = 1. (3.2) 

Proof 



From (|2.35| ), we have 

oo 

f n {x) = 2<J)- 2 (x) J <j) 2 (z)(w(z) - S n )f n ^(z)dz (3.3) 

x 

and , because of the equivalent form ( |2.36|) , 



X 

f n {x) = -2<p~ 2 (x) f 4> 2 (z)(w(z) - S n )f n ^(z)dz. (3.4) 







Assume that ( |3.2|) holds for f n (x) when n = m — 1. Let 

w(x) — E m = at X = x m . (3.5) 



Since w'(x) < and because of (|2.37| ), (|3.5| ) has one and only one solution. Furthermore, 



w(x) — £ m is negative for x > x m and positive for x < x m . It follows then f'^x) < for 

x > x m on account of (|3.3j), and for x < x m because of (|3.4|). When n = 0, /o(l) = 1 which 
satisfies (|3~7l|)-(|3~2"|). The theorem is proved by induction. 
Theorem 3.2 



En > £i (3.6) 



for all n > 1. 
Proof 
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From (|2.28| ), we have 



_ \w{x)f n _ l {x)} 

&n &1 — r n / M ^1 

[fn-l[X)\ 

[(w(x) - £i)/ n -i(x)] 

[fn-l(x)} 

Let x\ be the solution of 



(3.7) 



w(x 1 )-8 1 = 0. (3.8) 

By using 

[w(x)-S 1 } = 0, (3.9) 

we may write ( |3~7| ) as 

S n -S 1 = [(W(X) - S 1 )(fn-l(x) - f n -l{Xl))}/[f n -l{x)}. (3.10) 

Note that w(x) — Si and f n -i(x) — f n -i(xi) are both negative for x > X\, and both positive 
for x < X\. It follows then £ n — E\ is always positive. 
Theorem 3.3 



where for g sufficiently large 



< Si = 0{-=) « 1. (3.12) 



2. Introduce 



<7 2 ' 



I = 2 / <jr\y)dy J <p\z)w(z)dz. (3.13) 





and 



oo oo 

2/..\j„. /" /2/ 



J = 2j y <P\z)dz. (3.14) 

i y 



For g sufficiently large, 



J<0(^)«1, J<0(^)«1. 



0( — J<0( — (3.15) 
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The proofs of (|3.11| ) - ( [3.15| ) are given in Appendix A, together with the upper bounds of Ex, 
J and /. 

It is useful to define 



/ n (0), f n (l) and E n 

to be respectively the maxima of f m (0), / m (l) and £ m for all m < n. 
Theorem 3.4 



/n(0) < 



/n(l) < 



MX) 

i 

i - jz, 



and therefore 



En < 



E n < / n _ x (0) £ x < / n (0) Ex 



E x 1 E x 1 

< 



1 - J 1 - J£ n _i ~ 1 - J 1 - J£ n 



Proof From ( |2.36| ), we have 

f m (0) = / m (l) + 2 / <p-'\y)dy I cf> 2 (z)(w(z) - E m )f m „x(z)dz 



i y 

2/„,\j„, /" i2i 




1 3/ 



< / ro (l) + 2 / <T%)<fy / 4>\z)w{z)f m _x{z)dz 











On account of Q, /m-i(O) > / m -i(4 ( |3~2lD leads to 



/m(0) < frail) + /^(O) 2 j J f {z)w{z)dz 



/ m (l)+/ m _ 1 (0)/. 



(3.16) 



(3.17) 

(3.18) 
(3.19) 

(3.20) 



(3.21) 



(3.22) 



From the definition flHSft , / m (l) < / m -i(0) < / n _i(0) for all m < n. Thus ( |3~22| ) 

implies that for all m < n 



fm(0) < fn(l) + fn-l(0) I, 



and therefore 



/n(0)</„(l) + /„-l(0) /, 

16 



which, on account of / n _i(0) < / n (0), gives ( p,17|) 
From and ( CT) , we have 



For x > 1, 7/7(2;) < | 



£ n > St > - . 



< £ n - 7/7(2;) < £ n 



00 00 

2/„U., /" J.2/ 



and, on account of (|2.35| ) and ([ 

f n (x) < l + 2£ n I (f)-' 2 (y)dy I ^(z)f n _ 1 (z)dz, 



00 00 

-2/„,\j„, /" ^2/ 



< 1 + 2£ n / n _ 1 (x) / (f>-\y)dy / <^(z)dz. 



Set x = 1, we have 



/ n (l) < l + £ n / n _i(l)J 



which leads to 



/„(!) < 1 + En /n-l(l) J < 1 + £n f n (l) J 



and (3.18). 

Using ( gjp and (U), we obtain 

/n-i(0)H 



/n-l(O) 5; 



1) 



(3.23) 



(3.24) 



(3.25) 



/«-i(oo)[l] 

which gives (3.19). Combining (|3.17| ) - (3.19), we derive ( |3.20|) and complete the proof of 
Theorem 3.4. 

Let K be the lesser valued solution of the quadratic equation 

£1 1 



K 



l-Il-JK 



i.e., 



K = —(1-^1-AJ£ 1 /(1-I)). 



(3.26) 



(3.27) 
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For g sufficiently large and using ( |3.15| ), we find 



therefore 



^= T ^ 7 (i + I # 7 + ---); (3.2* 



K>£ 1 (3.29) 



and 



K — > E\ when g — > oo. (3.30) 



Theorem 3.5 



1. £ m < K for any m > 1, (3.31) 

2. lim £ m = exists and is < K. (3.32) 



Proof Assume (|3.31| ) holds for m = n — 1, we have from ( |3.20|) 

7- ^ 1 



From ( |3.29|) and the definition ( |3.16| ) of S n , we find 

£~i <~£~ 2 < ■ ■ ■ <T n < K 

for all n. Theorem 3.5 then follows. 
Theorem 3.6 For g sufficiently large, 

1. lim f n (x) exists for all x > 0. 

2. lim S n exists. 

n— >oo 

The proof is given in Appendix B. With Theorem 3.6, we establish the convergence of the new 
iterative solution. 

As we shall see, from a practical point of view this new iterative process is also quite 
effective in producing good approximate solutions of the quantum double-well problem. The 
first iteration fi(x) will be examined in Appendix D. 
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4. Generalization 



In this and the following sections we will extend our analysis to the lowest odd eigenfunction 
tp od . As an intermediate step we shall first construct a solution ip + of the Schroedinger equation 

(T + V — E + )ij + = (4.1) 

in the region 

x > 0, (4.2) 

where H = T + V is the same Hamiltonian given by ( |2.1|) , and ip + satisfies the boundary 
conditions: 

at x = 0, ^±( ) = ^±(0) = s-l (4.3) 

W+ <P+ 

and 

at x = oo, ^+(oo) = 0. (4.4) 
As we shall see, the solution 4>+(x) can be derived from 

f x ) = e -9So(x)^_ / 45 ) 
1 + X 

given by (2.6), through the Green's function 

D + = -2<f) + 6(f)- 2 6(f) + (4.6) 

and the integral equation 

il> + = <f> + + D + (u-e + )i{> + , (4.7) 
where u — (1 + x)~ 2 and 6 are given by ( |2.11|) and ( |2.16|) and f + is related to E + of ( ft.l|) by 



^+ = ^-^+, (4.8) 



analogous to ( |2.21|) . In this section we present the convergent iterative process leading from 



<f) + to ip + . In the next section we will show how ip + can in turn lead to the lowest odd 
eigenstate vp od of H. 
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As in ( gjp - we write 

^+ = 0+/+ , (4-9) 

where / + satisfies 

f + = l + D + (u-S + )f + (4.10) 

with 

D + = 0; 1 J D+0+ = -29(f)- 2 e(t) 2 + . (4.11) 
Differentiating (|4.9|), we have 

■0' 0' f' 



^+ 0+ /+ 

The boundary condition ( fOf ) implies that at x — 0, /+(0) = 0, and therefore 



oo 

| $.(«-£ + )/ + tte = 0; (4.13) 



i.e., 

oo 



/ u(f) + ip + dx 

= ^ • (4-14) 

o 

Analogous to ( |2.26|) , we introduce the iterative sequences {/+,«} and {£+,«}, with 

/ +i0 = 1 and £ +i0 = 

when n = 0. For n > 1, 

/+,„ = 1 + D+(it - e+, n )f+, n -i (4.15) 

and 

£+,n = [U / [/+,n-l] + , (4-16) 

where for any function F{x) 

oo 

[F] + = f <f)l(x)F(x)dx. (4.17) 



o 
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Because of (|4.16| ), 



4n(0) = 0. (4.18) 

Thus, ( |4.15| ) can be written either as 

oo oo 

U, n (x) = 1-2/ r+(y)dy J <f>l(z)(u(z) - £ + , n )f + , n ^(z)dz (4.19) 

x y 

or equivalently, 



x y 



/+,„(*) = /+, n (0) - 2 y ^(y)dj/ J r + (z)(u(z) - £ +>n )f +;n ^(z)dz (4.20) 



Parallel to the theorems given in Section 3, we establish the corresponding properties 
satisfied by and £ +>n in terms of Theorems 4.1-4.6 proven below: 
Theorem 4.1 At all x > 0, 



and therefore 



W) ^ ( 4 - 21 ) 



Z+,n(0) > f + ,n{x) > f +in (00) = 1. (4.22) 



Theorem 4.2 



S+,n > (4.23) 

for all n > 1. 

Proof By following ( |3~3"1 ) - ( |3 . 1 0| ) and changing /„, £ n , [F], w to S +t7l , [F]+ and u 
respectively, we can readily establish the above two theorems. Note that u' < 0, like w'. 
Theorem 4.3 

1. S +A = \ + ^ + ai7 (4.24) 



where 



311 , ln ^ 
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and, if we neglect 0(e sf), a\ is positive and 0(g 2 ) 
2. Introduce 



and 



We have 



i y 

2 J r+Wy J <p 2 + (z)u(z)dz (4.26) 





00 CO 

J+ = 2 j <j>+\y)dy J <j>%{z)dz. (4.27) 

1 y 



h < ~(ln M + 1) (4.28) 
g V 3 



and 



J+ < ^-ln (e + 2e^~g ). (4.29) 

^g 

Proof To prove this theorem, we need several formulas, already established in Appendix A 
for the proof of Theorem 3.3 in Section 3. Eq. (|4.24| ) is identical to (A. 2), and (A. 3) gives 

From (|4.26| ), we see that 



\2 



I + = 2jdyjdz e -MSo(z)-s m (1+^1 ? (4 30) 

I + z ) 

which is \ of the righthand side of (A. 60). By following (A. 61) - (A. 64), we see that 1+ is 
less than \ of the bound (A. 65), which leads to (|4~28| ). 

According to (|2.8| ), 4>{x) oc (f>+(x) for x > 1. Thus J+ above is the same J given by ( |3.14|) . 
Its upper bound ( |4.29|) is the same one given by (A. 39. 2). Theorem 4.3 is then proved. 

As in ( EnBD , define 



/+, n (0), /+, n (l) and £ + , n (4.31) 
to be respectively the maxima of / +>m (0), / +jm (l) and £ + , m for all m < n. 



22 



Theorem 4.4 



/ + ,„(0) < (4.32) 
1 



1 ~~ ^ + £ + ,n 



£+,n < / + ,„-l(0)5 + ,i < /+,n(0)£ +1 l (4.34) 



and consequently 



^<T^ ± 7 ^^^7^7 (4-35) 

1 — J + 1 — J + £ +;n _i 1 — i + 1 — J + C +)n 



Analogous to ( |3.26| ) - (|3.27|) , introduce 



K + = ^-(1 - y/l - 4J + 5 +il /(l - /+) ), (4.36) 



which is also the lesser valued solution of the quadratic equation 



Thus, K + > £ + i and — > when g — > oo. 
Theorem 4.5 



1. £+, m < fsT+ for all m > 1 (4.38) 



2. lim £+ m exists and is < K+ . (4.39) 

m— >oo ' 

Theorem 4.6 For g sufficiently large, 

1. lim f+, n (x) = f+(x) exists for all x > 0. (4.40) 

2. lim £ +n = £ + exists. (4.41) 

n— >oo 

The proofs of these theorems are identical to those of Theorems 3.4, 3.5 and 3.6, provided 
we replace f n , £ n , w, J, J and K by / +>n , £+,„, w, J + , J+ and K + . Thus, we establish 
the convergence of the new iterative solution for ip + . 

It is straightforward to verify that, similar to ( |4.21|) and (|4.22| ), for x > 



f+{x) < (4.42) 
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and therefore, 

/ + (0) > f+(x) > f+(oo) = 1. (4.43) 

Furthermore, 

/ + (0) < (1 - L)~\ (4.44) 
where L < In ^fg) is given by 

L = 1 - (1 -/+)(! - J+£+) = I+ + J+£+ - I+J+S+ (4.45) 



with bounds of I + , J + and £ + given by (|4.28| ), ( |4.29|) and S + < K + , on account of ( 



In addition, as mentioned before, after ( ft.!2| ), at x = 

/;(0) = 0. (4.46) 

Eq.(4.42) gives an upper bound for f' + (x). It is also useful to set a lower bound for f' + (x). 
Differentiating (4.10), we have 

oo 

-2/„\ [ ±2 



f + {x) = 2<p-\x) J r + (z)(u(z) - £+)U(z)dz. (4.47) 

x 

Because of (4.13), the same expression can also be written as 

X 

f'Jx) = -20; 2 (x) f <P 2 + (z)(u(z) - £+)f + (z)dz. (4.48) 





2 



Since u{z) = (1 + z)~ 2 > 0, (4.47) gives 



oo 



f' + (x) > -2S + ^-\x) / 4>\{z)U{z)dz 



> -2£ + U(x)<p-\x) J 4>l(z)dz (4.49) 

X 

From (4.23)-(4.25),(A.4)-(A.5) in Appendix A and neglecting 0(e~l 9 ) we have S + > \ . The 
replacement of £ + by \ in (4.48) yields the inequality 

x 1 
f + {x) > -2d>l 2 (x) J ' <f + {z){u{z)--)U{z)dz 



x 1 
> -2/ + (O)0; 2 (x) I <f>l(z)(u(z) - -)dz. (4.50) 
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For x > 1, 05- (x) = e" 2950 ^^^) 2 < e' 29So( - x \ Also, according to (A.33)-(A.36) proven 
in Appendix A, 

00 

j(x) = e 29So ^ [ e- 29S °( z) dz < — — — — — (4.51) 
J (1 + x) 

where 



C 



+ -£= + •••) ■ (4.52) 
V 9 V 2^ ^ 1 J 



Consequently, for x > 1, (4.49) leads to 



f + (x) > ~£ + f + (0)C . (4.53) 



For < 2; < x < 1, we have 

So(z)-S (x) = (x-z)(l-\(x 2 + xz + z 2 )) 
> (x - z)(x - -(x + z + z)) 



and 



which yield 



07 2 (x)0 2 (z)(u(z) --) = e 2fl(So(x)-SoW)(i±5 )2( _ 1) < 3e -ffl(— ) 2 . (4.54) 

+ v 7 +v /v v 7 4 7 1 + z v (l + ^) 2 4 7 v ' 

Thus, from (4.50) and for x < 1, we derive 

X r— 

> -6/ + (0) / e-l^-) 2 ^ > ~J—U(0) ■ (4.55) 



Combining (4.52)-(4.53) with (4.55) and for g sufficiently large we have, at all x, 

/+» > —o\h f +^ ^ 
z V 9 

which is a convenient (though not the best) lower bound. 

In Appendix D, we will give an analysis of the first iterative solution / +il (x). 
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5. The Lowest Odd Wave Function 



5.1. Formulation 

The lowest odd wave function tp od of the double-well potential V — ^-(x 2 — l) 2 satisfies 
a similar Schroedinger equation as ([2.11): 

(T + V — E od )ip od = 0. (5.1) 

with 

^od(-x) = -ipod{x)i (5.2) 

and therefore, ^ o d(0) = 0; besides at x — 0, ip od has no other node. In this section, we give 
the derivation of ip d(x) through an alternative convergent iterative procedure which differs 
from those of the previous sections. In this new approach, as we shall see, each successive 
iteration yields an improvement by a factor 



4 



e = e~-- 9 (5.3) 

that is much smaller than the previous 0(g~ l ) factor. 

We begin with the function ip+(x) of (|4.1|), that has been derived in the last section. 
Recall that while 4>+(x) satisfies the differential equation (T + V — E + )i}) + = 0, it is not an 

eigenstate of the Hamiltonian H = T + V , since at x = —oo, ^+ diverges. However, as 
will be shown below, it is possible to construct a good trial function x{ x )> °dd in x and well 
behaved everywhere, including x = ±oo. For x > 0, we write 

X(x) = ip + (x) - i/j-(x), (5.4) 

where 

oo oo 

il>-{x) = 2^ + {x) J ^ 2 {y)dy J ^ 2 + {z)dz (5.5) 

x y 

with the constant 7 chosen to make \ vanish at x = 0, 

x(0)=V+(0)-V>-(0) = 0. (5.6) 
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For x < 0, we set 

X(x) = -x(-ac) (5.7) 

Because of (|5.6| ), 7 is a positive constant given by 

00 00 

r l = 2jr + \y)dyj^l(z)dz. (5.8) 

y 

It can be readily verified that at x = 0, both and x'(a;) are continuous, but not 



x"( x )- Furthermore, by differentiating (|5.5| ), we have 

(T + V- E + )ijj„ = -7^+. (5.9) 

Define u(x) by 

v(x)x{x) = -<y4> + (x) (5.10) 

for x > 0, and u{x) = v(—x) for x < 0. 
Theorem 5.1 

(i) x is the lowest odd eigenstate of the Hamiltonian T + V + v, with eigenvalue E + ; i.e., 

(T + V + z/-£ + ) X = 0. (5.11) 

(ii) 

u(x)<0, (5.12) 

and for x > 

> 0. (5.13) 

(Since z/(x) is even in x, v'(x) < for x negative.) 
Proof 



(i) Write as 



X = ^+(1-t z )- (5.14) 
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From 



C^-)' < 0- (5.15) 
Since at x = 0, ip_ = we have for x > 

V>-(a:) < V+fc), (5.16) 

and therefore 

X (x) > 0, (5.17) 



in which x( x ) = occurs only at x = (and x — > oo). In ( |5.15|) -( |5.17|) , x is restricted to 



x > 0. (In the following, unless specified, the same restriction applies.) 

From ( |5.7|) and ( |5.9|) -( j5TT0| ), we see that xi x ) satisfies the Schroedinger equation ( |5.11|) 
at all x. Since x( x ) nas or| ly one node at x = 0, it is the lowest odd eigenstate of the 
Hamiltonian T + V + v. 

(ii) From ( |5.10| ), we have v < 0. 

Multiply ( ft~H ) and ( |5.11| ) on the left by % and ip + respectively, and take their difference. 
We obtain 

Consequently 

'^+0*0 y _ . -2, v / , 2 



( W } = ~ 27x (x) / ^ y)dy - (5 - l8) 

Since the left side is —i/fj, we derive ( |5.13| ) and complete the proof for Theorem 5.1. 
Theorem 5.2 The constant 7 satisfies 

7 < AgJ— e"t 9 (l + a 7 ) (5.19) 

V 7T 



and 



7 >4W^e-t ff (l-/3 7 ), (5.20) 



28 



where 



a 7 = 0(-l nv ^) and (3, = 0{— -), (5.21) 

9 9 lhi 

so that for g sufficiently large 

7 4c/W— e"t 9 . (5.22) 

V 7T 

The explicit form and the derivation of these bounds « 7 and /3 7 will be given in Appendix 
C. Here we only discuss how to obtain the approximate formula ( |5.22| ). 
In (|5^), because of ( |Q3"D -( |3~33D , for g sufficiently large 



Hence, ( p.8| ) can be approximated by 

oo oo 

r x = i\ r+(y)dy J 4> 2 + (z)dz. (5.23) 
o y 

According to (2.6), the ^-integration is dominated by the region z — 1 + 0(-^=), where 

<Pl(z) = e~ 29 ^ 1)2 (5.24) 

is near its peak, whereas the y-integration is largely determined by the region y near 0, and 
O(i), where 



= (\±ly e by 3 -igy+h g* l e -*9v+ig. ( 5 . 25 ) 



Thus, we may futher approximate ( |5.23|) by 

oo cx_ 



' ' ; 4<? V2g 



-oo 

which gives ( |5.22|) . Therefore, on account of fl573]), 7 = 0(e), and so is v(x). 

Because x ar| d are both the lowest odd eigenstate of two closely related Hamiltonians 
which differ only by 

i/(x)=0(e), (5.26) 
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the iterative approach for the derivation of ip od by using \ as the zero' h approximation should 
be an effective one, as will be shown. (In the derivation, v{x) appears only as part of the 
product ux', therefore its effect is always small, 0(e), even though u(x) by itself is singular at 

x = 0.) 

The Green's function associated with X ' s 

G = -2 X 9 X - 2 9X, (5.27) 
where, as before, 9 is given by ( |2.16|) . Because of ( |5.11|) , G satisfies 



(T + V + i/-E+)G = l. (5.28) 
Therefore, the eigenstate ip od , defined by ( [STID , is the solution of 

*Pod = X + G{v + A od )ip od , (5.29) 

where 

Aod = E od - E +] (5.30) 



i.e. 



oo oo 
-2 



ipod(x) = x(x) - 2x(x) / X (y)dy / (v{z) + A od )x(z)ijj od (z)dz. (5.31) 



y 



At x = 0, in order that the above equation agrees with vp od (0) = x(0) = 0, the integrand 
must satisfy 



limx 2 {y) (v(z) + A od )x(z)ij od (z)dz = finite, 
y^o J 
y 

which implies 

oo 

J (u(z) + A od ) x (z)i> od (z)dz = (5.32) 
o 

and, on account of (|5.10| ), 

oo 

f ip + (x)ip od (x)dx 

Aod = 7^ >0. (5.33) 

/ x(x)*Pod(x)dx 
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The positivity of A od = E od — E + is to be expected. For x > 0, let ip be the solution of 
the second order differential equation 

(T + V - E)ip = 0. (5.34) 

Assume that tp has no node when x > 0, and it satisfies the boundary conditions: 

— = p>0 at i = 

i> 

and (5.35) 
■0 = at x = oo. 

When the slope p = 0, ip = ip ev and E = E ev . As p increases, so does E. When p = 

(j>' + (0)/4>+(0), tp — and E = E + > E ev , and when p = oo, ip = ip d and E = E d, with 

E od = E + + A od > E+. (5.36) 

Write 

ipod(x) = x(x)k(x). (5.37) 
On account of (|5.31| ), k(x) satisfies 

oo oo 

k(x) = 1 - 2 J X - 2 (y)dy J Hz) + A od ) X 2 (z)k(z)dz. (5.38) 

x y 

Define 

G = -26 X - 2 e X 2 . (5.39) 

Eq.( |5.38|) can also be written as 

k = 1 + G(v + A od )k. (5.40) 



As in Q2.26|) -( [2~2"8"D and ( |4.15| )- (|4.17| ), k(x), and therefore ip od (x), will be solved by an iterative 
procedure. 

Introduce the sequences {k n } and {A n }, with 

k = 1 and A = 0. (5.41) 
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For n > 1, 

fc n = l + G(z/ + A n )fc n _! (5.42) 

and 

A n = -{z/^x}/!^-!}, (5.43) 

where for any function F(x) 

oo 

{F} = J X 2 (x)F(x)dx. (5.44) 



o 



Thus, ( |5.42| ) can be written as either 



oo oo 

k n (x) = 1-2 [ X ~ 2 (y)dy f Hz) + A n ) X 2 (z)k n . 1 (z)dz } (5.45) 



x y 

or 



x y 



k n (x) = fc n (0) - 2 y _/ M*) + A n ) X 2 (z)k n ^(z)dz. (5.46) 



5.2 Convergence 
Theorem 5.3 

(i) For n > and x > 0, we have 

> 0, (5.47) 

and therefore 

K{0) < k n {x) < k n (oo) = 1. (5.48) 

Furthermore, 

k n {0) > K, > (5.49) 

where K, is independent of n. 

(ii) For n > 1, 

Ai > A n > 0. (5.50) 
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Proof Assume that (|5.47| )- (|5750| ) hold for n = m - 1. From ( |5.45|) and ( |5.46|) , we see that 

for n = m, 

oo 

Kn( x ) = 2 X' 2 (x) J(u(z) + A m )x 2 (z)k m _ 1 (z)dz 

X 

X 

= -2 X - 2 (x) J(v(z) + A m ) X 2 (z)k m -i(z)dz. (5.51) 



By assumption, k m -i(x) > 0. From (|5.12| ), v < 0; therefore, ( |5.43|) gives A m > 0. 
Let 

u(x) + A m = at x = x m . (5.52) 

Since, according to ( |5.13|) , v' > 0, we see that ( |5.52| ) has only one solution. Furthermore, 
u(x) + A m is positive for x > x m , and negative for x < x m . Thus, by using either the upper, 
or the lower equation in ( |5.51] ), we find k' m (x) > when x > 0. 
From 

- A m + Ai = 1 v ; u n + Ai 
{fc m _i(a;)} 

{(u(x) + Ai)fc m _i(x)} 

{/c m _i(a;)} 

According to ( |5.52|) , z/(x) + Ai = at x = x\. Since k = 1, we have {u(x) + Ai} = 0. 
Write (|533| ) as 

{(z/(x) + Ai)(A; m _i(x) - fcm_i(a;i))} 



A m + A x 



(5.54) 



{fc m _i(x)} 

When x > x\, both z/(x) + A x and k m ^i(x) — fe m _i(xi) are positive; when x < xi, both 
negative. Therefore, for m > 1, A m < A x . To complete the induction, we still have to prove 
that when n = m, (|5~49| ) is true, also ( p~47|) -( |5~49l) hold for n = and ( |530| ) is valid for 

n = 2. 

By assumption, < k m -i(x) < 1. Since z/(x) is negative, and for m > 1, A m is positive 
but < A x , we have 

oo oo 

+ A m )x 2 (^)£: m -i(2)Gb < A m x 2 {z)dz 



< Ax / x 2 (^)^; 



(5.55) 



33 



therefore at x = 1, ( |5.45| ) implies 

00 oo 

k m (l) > 1 - 2A 1 J X ~ 2 (y)dy J X 2 (z)dz. (5.56) 

1 y 

From ( |5.46| ), we have 

i y 

MO) = M 1 ) + 2/ x~\y)dy J iy{z) + A m )x 2 (^ m -i(20^ 


1 2/ 

> M 1 ) + 2 J x~ 2 (y)dy J is(z) X 2 (z)dz, 



which, on account of vx = —l4>+, gives 



1 

2/ 



MO) > Ml) " 27/ X~ 2 (y)dy J X (z)^ + (z)dz. (5.57) 




Define 



and 



1 y 

1 = 2 J X ~ 2 (y)dy J x(z)Mz)dz (5.58) 


00 00 

J + = 2 1 x" 2 ^ I X 2 (^)^- (5.59) 

1 2/ 

We see that from (|jg)-(g5^), 

fc m (0) > 1 -7fj + J7 + — 1 e= £. (5.60) 

7 J 

In Appendix C, an upper bound will be given for J+ l 7 +: ^ L , which is 0(e°). (See (C.86), 
(C.87), (C.107) and (C.118).) Since 7 satisfies (pJ|)- (^22l) , we have /C = 1 - 0(e) > 0. 



When 71 = 0, k = 1 which satisfies (|5~47lH|5~49l ) . This leads to ( ^~47|) -( ^~49|) being valid also 
for n = 1; in addition, A 2 > Ai (since for m = 2, the proof (|5.53| )- (|5.54| ) depends only on 
ki satisfying (|5.47| )- (|5.49| )). Theorem 5.3 is then established. 



Theorem 5.4 



1. lim k n (x) = k(x) exists for all x (5.61) 

n^oo 

2. lim A n = A od exists (5.62) 



n— >oo 
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Proof The proof is parallel to that of Theorem 3.6, or Theorem 4.6. 
Introduce 



Kj(x) = k n+j (x) - k n (x) (5.63) 

and 

k ™ = max overall x. (5.64) 

As we shall see, 
Lemma 

K y-R Y 1 < r, (5.65) 

where r = 0(e) is finite and independent of n and j. 

In the following, from ( |5.66| )- (|5.91| ) we first give the proof of the lemma. 
According to (|5.42| ) and ( |5 .63|) , we have 



K 



3 



G((u + A n+j )k n+j ^ - (u + A n )k n A =G(Q 1 + Q 2 + Q 3 ) } (5.66) 



where 



Qi = (A n+j - A n )k n+j -i (5.67) 

Q 2 = A n (/c n+ j-i - fc n _i) (5.68) 

Q 3 = v{k n +j_ x - k n -i). (5.69) 

Define for x > 1, 

oo oo 

q a (x) = GQ a = -2 / X ' 2 (y)dy I X 2 (z)Qa(z)dz (5.70) 

and for < x < 1, 

q a (x) ee 9a (i) + 2 y y X 2 (z)Q a (z)dz, (0.7] ) 

x o 

where the subscript a = 1, 2, 3. Because of ( |5.45| )- ( |5~46f) and (|5.66| ), we see that 



\ v 

2/„ \ t . [ 2/ 



k?(x) = gi(x) + g 2 (x) + g 3 (^) (5.72) 
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for all x. 

On account of ( ggp and (|5T63| ) , 



A --A = i^Zll Mfcn-l + Kj X )} 



J- 

{z/fc n _l}{fc n _l + tf- 1 } - {Vjkn-l + 4~ 1 )}{fen-l} 

{k n +j-i}{k n -i} 



f] , r, , (5-73) 
From (|5.48| )- (|5.49| ), we have k n+ j-i(x) > /C; consequently, the absolute value of the above 



expression satisfies 

By using the definitions (|5.64| ) for k ™ and ( |5.43| ) for Ai (when n = 1), together with 
( |5.50| ), we derive 

IWJ- 

{1} 

K]-\A n + A 1 )/)C 



\A n+j -A n \ < «:^ 1 (A rt + ^)//C 



< 2A 1 K r ;~ 1 K . (5.74) 
Since < k n+ j-i(x) < 1, from (|5.67| ), it follows that 

|Qi| <2A 1 k]- 1 /)C, (5.75) 

and, because of ( |5.70| )-( p71] ), 

\q 1 {x)\<2JA l K n f 1 /K (5.76) 

where 

J = J + + J- (5.77) 

with J+ defined by ( g3g ) and 

i y 

J' = 2 J X~\y)dy J x\z)dz . (5.78) 
o o 
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Next, we examine Q2 given by (5.68). Note that < A n < A\, according to ( |5.50|) 
Therefore, by using (|5.63| )- (|5.64| ), we fined 



\Q 2 \ < A{K n f\ 



(5.79) 



which, on account of ( |5.70|) -( [5~7"l~D , leads to 



\q 2 (x)\ <JA 1 k]~ 1 



(5.80) 



For Q 3 , multiply its definition (5.69) by X- Using (|5.10| ), we have 



XQ3 = - kn-i) . 



(5.81) 



Both 7 and ip + are positive. Thus, 



\xQs\ < 7^+ K 



n-1 



(5.82) 



Combining this expression with ( |5.7C1| )-( |571] ), we find for x > 1 



00 00 

..ts n-1 f ,,-2/ 



\q z (x)\ < 2^2~ J X~\y)dy J X (z)^ + (z)dz 

x y 



and for < x < 1 



(5.83) 



1 y 

\q 3 (x)\ < |g 3 (l)| +2 7 ^r 1 / X~\y)dy J X (z)^ + (z)dz 

x 

As will be established in (C.53) of Appendix C, for x > 1 



(5.84) 



ip-(x) < 7J+^+(x) 



(5.85) 



where J + = 0(-ln^/g) is bounded by ( |4.29|) ; therefore 



x ( x ) = ip + (x) - ij}_(x) > (1 - -fJ+)ip + {x); 



i.e. 



x -\x)<{\- 1 j + r^-\x). 



(5.86) 
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This together with x( x ) > 4'+( x ) g' ve - f° r x > 1, 

oo oo 

\q 3 (x)\ < 2 7 «f x (l -lJ + y 2 J r + 2 (y)dy f ^\{z)dz . (5.87) 

1 

Note that for y < 2 



Since, according to (|4.27| ), 2 / <j)^ 2 (y)dy f 4> 2 + {z)dz = J+, we derive for x > 1 

1 2/ 

|g 3 (x)| < 7 J + (l-7^r- ( 5i 
From (|5.84| ) and the definition of X given by ( |5.58| ), it follows that, for x < 1, 



i ?3 (x)i <7 [x+ j+(i- 7 ^+r 2 ] «r' • ( 5 - 89 ) 



Combining together (|577g| ),( F^DD , (p^5|) and ( |CT| ), we obtain 



where 



<lA = r (5.90) 

A = J + J + (l - 7 J+r 2 + J— (1 + p), (5.91) 
confirming ( |5 .65|) . The bound of J + = oQlny^) is given by ( ft.29| ). Explicit bounds for 

l = 0(g-l), J = 0(-hi^g)+0(g-l) , 

^1^1 and j9i4gj— e '^ 9 (5.92) 
7 V 9 

are given by (C.86)-(C87), (C.109), (C.118) and (C.1)-(C8) of Appendix C. The bound for 
K, is given by ( |5.60|) . The lemma is then established. 
From (|5.48j ), we see that when n = 1, k j satisfies 

k j = max |/c 1+ j(x) — ki(x)\ < 2 . (5.93) 

Thus, (ET90D gives 

K] < rTc n f l < r 2 n ]~ 2 < ■■■ < r n ~ l K ) < 2^ . (5.94) 

Analogous to (B.24WB.25), lim A; n (x) exists for all x and likewise lim A n exists. Theorem 

v ' v ' n— >oo v ' n— >oo 

5.4 is proved. 
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Appendix A 



In this Appendix we give the proof of Theorem 3.3 of Section 3, together with some related 
inequalities. The bounds for £ 1( J and / of (|3.11| )- ( [3.1 5|) will be given below in Sections (i), 
(ii) and (iii). 

(i) We shall first establish the following three statements (A.l), (A.2)-(A.3) and (A.4)-(A.5). 

oo 

J <f) 2 {x)g(x)dx < 2(2g + 3)e~^. (A.l) 



Write 



/ 4> 2 + (x)u(x)dx 

^ = z + 56 - + ( A2 ) 



°1 < 7^2- ( A3 ) 





then 

311 

2V 

Furthermore, if we neglect 0(e~^ 9 ), a\ is positive. According to ( 3.11|) , the first iterated 
energy E\ can be written as 

ft-i + ii + ft. (A4) 

As we shall see, 

= ai + 0(e-3 s ). (A5) 

Proof From (2.6) and (2.7), we have 



(A6) 



For < x < 1, because of and (gl5|) , 



and 



therefore, 



5-+ 10+/ 5- + 1 0+ 



(i + (ttt)F); ( a ^ 



{l + xf v v #+l 



^ = TTtS^ 6 "^ 1 + (^^"H (A9) 
(1 + xy g + 1 g + 1 
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Since g = for x > 1, the left side of (A.l) equals 



i 



(x) 2 g(x)dx = Ag^— -)e-i 9 (l + r), 
9+1 



where 



On account of (|2lp , 



r = 2 (^4) [^—e^-l°dx. 
g+l J (1 + x) 2 



For < x < 1, we have 
and 

1 1 

4 < (1 + z) 2 < ' 



Hence, 

1(1 - e^) < f j-^e^^dx < A(i _ e -|i 
8# J (1 + x) 2 4g 

Combining (A.lO)-(A.ll) and (A. 15), we derive 

9~ 1 ^ -tar, , 3 ,9- 1 



/ 2 (x)^(x)rfx < 4^(^— -)e-^{l + — -)(1 - e 
7 # + 1 2c? g + 1 



and 



cxj 

U\x)g{x)dx > 4g(i^-) e -i'{l + U 9 —^)(l ~ e 

J + 1 4(7 + 1 



,0—1, 4„ r 1.0—1 

54— e-a^{l + r ^— r 
+ 1 ^g g + l 

u 

The inequality (A. 16) implies (A.l). 

Next we turn to (A. 2). From ( [2.11] ), we have 

1 

.2(1 + x) 2 ' 2 2 (l + x) 
By partial integration and using (2.6), we obtain 



U(x) = — (X — 1) ( — r— + 

V ; 4 v J \2(l + x) 2 



(f) 2 Jx){u(x) - -)dx = -e~ 2g5 °l _, M s M x ' 



4 y ' 1 2^(1 + a;) 2 v v ' 4 

+ / e - 2 ' S °{2^(T^M*>-I>}' 
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which, on account of (A. 18) and S ; Q = x 2 — 1, leads to 

oo oo 

/£(,)(«(«) - \)d* = + -/e^.(_ + (A20) 

Next, write 

oo oo 

/ 02 9 d = 1 /■ 2g s 1 / 5 2x 

v ' 

Taking the difference (A. 20) minus (A. 21) and doing another partial integration, we derive 
(A.2), in which 

1 f 1 , 103 1. 



+ i e V (1 + a;) 5 I (1 + a;) 3 + 2(1 + xf + 2 2 (1 + x) + 2V J' ( } 
in which the integral in the numerator is less than 

7 2o <? / 2 N2 1 , np . 54 45 36, , 311 1 f , 2 , , , , . n . 







Therefore, (A. 3) follows. Furthermore, if we neglect the 0(e 39) term in (A. 22), a% is positive. 
To establish (A.4)-(A.5), we define 



00 00 

2j„ f a2 



A = I <p z dx - / <p\dx (A24) 

and 



00 00 



B = J (f) 2 udx - J (j) 2 + udx. (A25) 


From (P78|), one can readily verify that 

00 

A = 2( ^)e-f»{ (1 + i(|^)e-|»)/^ 

1 



+ 2 + 2 (^)/(T^^-^} 
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and 



-1,1-1 00 



"9 + 1 
+ - + 2(- y 



udx 



6 9+1 { (l + x) 4 
Since for x > 1, < e -2 ^ -1 ) 2 , we derive 



i nr 



By using (A.13)-(A.15), we have 



and 



udx < 



1 FT 
8V2? 



/ —±—e 2 ° So -l°dx < j —±—e 2 ° s °-l°dx < A 
7 (l + x) 4 7 (l + x) 2 4o 



It follows then 



4( 



g ~ 1 
9 + 1 



e-* 9 < A < Ae~* 9 \ 1 



i rv 



(1 



4g 



and 



Therefore, 



Since 



3 V ^ + 1 



)e _ 5» < B < -e"3 9 U 



3 /~f 

28V 2^ 



(l + - e -3fl) + 



7^7 



A = 0(e~39) and B = 0(e~3 s ) 



/ (j) 2 (u + g)dx 
c _ o 

£ 1 — 55 

/ <p 2 dx 





we have 

oo oo 

b\udx + B + J <p 2 gdx 

c _ 

C l — oo ) 

/ <\)\dx + A 
o 

From (A.l), (A.2), (A.27) and (A.28), we see that (A.5) holds; i.e., 

= ! + ! I + ai + 0(e-S'). 
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with ax < 311/2V. Thus we establish fglTD -( gl2D of Theorem 3.3 



(ii) To set a bound on the integral J of (|3.14| ), we first introduce 

oo 

j{x) = e 29So{x) J e- 2gSo(z) dz (A32) 

X 

and then show, for x > 1, 
where 
and 



j{x) < t-^, (A33) 



(l + x) 2 

1 l + l 
2~g T~l 



C=— l -±± (A34) 



Thus, 



and, (A. 33) gives 



2/ng. (A35) 



C=,Rl + -i = +...) (A36) 



9 2v/2^ 



j(x)dx < -C. (A37) 
i 

As we shall see, the same C is also an upper bound for J: 

oo oo 

J = 2 [ <f)- 2 (y)dy f <j) 2 (z)dz < C. (A38) 



Next, we shall improve this bound by proving 

J < <'-V5 + £' n £ (A39.1) 

1 



< —l n (l + 2y/2ng) + —. (A39.2) 
where / is the same one given above. 

Proof For x > 1, S^x) = |(x - l) 2 (a; + 2) > (x - l) 2 and therefore at x — 1 



OO 

j(l) < / e" 2 ^- 1 ) 2 ^ = -y^. (A40) 
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Differentiating (A. 32), we find the derivative of j(x) to be 

f(x) = 2g(x 2 -l)j(x)-l. (A41) 
Using (A. 35) to define I, we see that I satisfies 

2^-1)1^-1 = 0. (A42) 
It follows then, for 1 < x < I, j'(x) < and therefore 

fa)<ft)<\^Yg (A43) 
Furthermore, partial integrating (A. 32), we find 

oo 

which gives, for x > 1, 

>'<*> < ^3T)' < A44 > 

The constant C in (A. 33) is determined by setting 

C 1 



(A45) 



( Z + l)2 2^-1)" 
For x > I > 1, since 

/ + 1 x + 1 

n^r > x- 1 ; 

we have x 2 - 1 > (1 + a;) 2 (7 - 1)/(Z + 1), and because of (A.44)-(A.45), we obtain 

< 2,(1 + ^-1) = (lf^' (A46) 
which becomes, at rr = /, 

^ < (TTiF = WE ' (A47) 

on account of (A.42) and (A. 45). Combining (A.43) with (A.46)-(A.47), we prove (A. 33). 

For z > y > 1, 

0- 2 (y)0 2 ( 2 ) = e 2 5 5o( J /)-2 s 5 (,) ( i±^)2 

\. ~\~ Z 

< e 2gS (y)~2gS (z) _ (A48) 
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Thus, (A.37) implies (A.38). 

To improve the bound (A. 36), we observe that 

S (z)-S (y) = ( Z -y)^( Z 2 + Z y + y 2 )-l) 

= (z-y)(^(z-y) 2 + y(z-y) + y 2 - 1). 

Thus, for z > y > 1, 

S (z)-S (y)>(z-y) 2 (A49) 

and also 

S (z)-S (y)>(z- y )( y 2 -l). (A50) 

Next, we separate the integration range in (A.38) for J into two regions: 1. 1 < y < a, and 
2. y > a; correspondingly, 

J = Ji + J 2 , (A51) 

where 

a oo 



J l= 2 j dy j dzr\y)<p\z) 



i y 
and 

oo oo 

J 2 = 2jdyJ dz<f)- 2 {y)(f) 2 {z). (A52) 

a y 

From (A.48)-(A.49), we have 



a oo 



-h<2 j dy j e~ 29(z - y)2 dz = (a - 1)^ (A53) 

and 



1 2/ 



oo oo oo 



J 2 < 2 (dy ( e-W*-M*-*dz = -f -P- = -ln^. (A54) 
J J q J y 2 — 1 2q a — 1 

Thus, for any a > 1, 



q J y 2 — 1 2o a — 1 



j <( tt - i Vl + ^ in ^' (A55) 

which is minimum when 

a 2 - 1 = ./A. (A56) 
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Comparing the above equation with (A. 42), we see that 



a = I 



\ 



2 1/2 
1 + J— = 1 + -J— + 
V irg 2 V ng 



(A57) 



consequently, (A. 39.1) is established. By setting a = 1 + (2ng)~^ , instead of/, we obtain the 

second bound (A. 39. 2) for J. 

(iii) For any function F(x), define 



i v 

2/.\i. f ±2, 



1(F) = 2 / r\y)dy / <j)\z)F(z)dz. 



(A5£ 



The integral / of (|3.13| ) can be written as 



/ = I( W ) = I( U ) + 1(g) 



(A59) 



with w, u and g given by ( |2.14j ), (|2.11j ) and ( |2.15|) . For I(u), we have, on account of ( |2~8"1 



i y 



I(u) < 8 / dy / ^ e - 29{SoW - So(y)) 



(l + y) s 
(! + *)< 



(A60) 



o o 

Next, divide the integration range into two regions: 1. 0<z<y<b<l and 2. 
b < y < 1 and < 2 < y. In region 1 

Sb(«) - S (y) = (y-z)(l- \(y 2 + yz + z 2 )) > (y - z)(l - ±(y + y + 1)) (AQ1) 



and in region 2 
In both regions 



SoO) - S (y) > (y - z)(y - ^(y + z + z)). 



(1 + ^) 4 



(A62) 
(A63) 
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Thus, (A.60) leads to 



b y 

I(u) < 32 J dy J dz e -!fl(»-*)( 1 -») 


1 V 



+32 J dy J dz e -l 9iy ~ z)2 



b 



32 / dy 



(l_ e -|ro(i-i/)) 



1 2/ 



+32 J dy J dz e -l 9iy - z)2 
b o 

< ffl{-^'»(i-») + i(i-»)^f 



whose minimum is at 



we obtain 




6 = 1 



/W<-(i"/f+i). 

V 3 



(A64) 



(A65) 



To analyze 1(g), we start from (A. 9). Since (1 + x) 2 < 1, (g — l)/(g + 1) < 1 and for 

< x < 1, e 2 9 s o-^g < i we have 



From (2.6)-0, we have 



2 (x)(?(x) < 16g e * 9 . 



(A66) 



which, on account of (A. 13), is < e * 9X for < x < 1. Thus 



I(# = 2 / 0- 2 (2/)A/ / 2 (z)£(^ 



/" 18 
< 32g / e-l 9y ydy < —. 
J Q 



Combining (A. 65) and (A. 67) we derive 



TV \ 24. , 7. 
/ = /(«,)< —(in Yy + j)- 



(A67) 



(A68) 
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Appendix B 



To establish Theorem 3.6, we introduce 

5]{x) = f n+j {x)-f n {x) (B.l) 

and 

A™ = maximum of \6j(x)\ overall x > 0. (B.2) 

As we shall see, 

Lemma For g sufficiently large, 

■Ml<R«l, (B.3) 

j 

where R is independent of n and j. 
Proof From ( |2.27| ) we can write 

Sj = D((w - £ n+j )f n+j -i - (w - £ n )fn-l) 

= D(L 1 + L 2 + L 3 ), (BA) 

where 

L\ = (£ n — £n+j) fn+j-i 

L 2 = (\-£n)(fn +J -l-fn-l) (B .6) 

Ls = (w-^){f n+j ^-f n ^). (B.7) 
Define for x > 1 

oo oo 

X a (x) = DL a = -2 J r\y)dy J <p 2 (z)L a (z)dz (5.8) 

x y 

and for < x < 1 

i y 

2/..U.. f a2/ 



X a (x) = A (l) + 2 J <tr\y)dy J <j) 2 (z)L a (z)dz } (B.9) 

x 

where the subscript a — 1, 2, 3. Because of ( |2.35| )-( |2~3~7l) and (B.4), we have 



5]{x) = X 1 (x) + \ 2 {x) + \ 3 {x) 
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(5.10) 



for all x. 

To set a bound for \\{x), we observe that, on account of ( |2.28|) and (B.l), 



p _ p [w/n-l] Hfn-1 + 



[fn-l] [fn+j-l] 
w/n-l][/»-l + - [wUn-l + K^Wn-l] 



[fn-l)[fn 



+3-11 



[fn+j-l] ° [fn+j-l 

From ( |3.2|) , we have [/ n+3 -_x] > [1] and therefore the absolute value of (B.ll) satisfies 

\£n-£n + j\<^- ] (£n\[Sr 1 }\ + IH^W) 

By using (B.2) and (3.31) we find 

n-l 

\£n ~ £n+j\ < A" (£n + 



= A]- 1 (£ n + S 1 )<A]- 1 2K, (B.12) 

where K is given by ( ft.27|) . Since 1 < f n +j-i(x) < f n +j-i(0), from (B.5) it follows that 

|Li| < [S n -E n+j \ f n+j -i(0). (5.13) 
On account of (B.8)-(B.9) and (B.12), for x > 1 

/OO y CO 

r\y)dy j <P\z)dz = 2 / n+i _!(0) K A]- 1 J, 

(5.14) 

where according to (A. 39. 2), 



for x < 1 

|Ai(x)| < |A 1 (l)| + 2/ n+3 -_ 1 (0)^- 1 2 j\- 2 {y)dy jy\z)dz 

= \\ x {\)\ + 2 / n+i _!(0) K A J "- 1 /(1), (5.15) 
where, in accordance with (A. 58), 

i(i) = 2 f l r 2 (y)dy f \\z)dz 

Jo Jo 
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<8 f'dy [" 'dz e -MSo(z)-s oi y))^ i 



J o * 7o (1 + z) 2 

<32 f dz e -MSo(z)-So(y)) < 24^/ M+i), (jB .i 6 ) 

Jo Jo g V V 3 y 

which is, on account of (A.61)-(A.64), the same upper bound for I{u) given by (A. 65). 
Combining (B.14) with (B.15), we see that at all x, 

\X x (x)\ < 2 / n+i _i(0) K(J + I(l))A]-\ (5.17) 

Next, we examine L 2 and L 3 of (B.6)-(B.7). Because of (3.6), (A.4) and (3.31), \ < 

S n <K, 

\i*\ < (£>--) < (k-±)aj-\ 

Thus, using (B.8)-(B.9) and by following the same steps leading from (B.13) to (B.17), we 
obtain 

|A 2 (x)|< (K-\)(J + I(l))A]-\ (£.18) 

Likewise, 

\L 3 (x)\<\w-\\A^\ 
Since, from C^TTTQ and ( ^TT4|) -( ^TT5l) , for x > 1, w(x) = u(x) < \ and for < x < 1, 



w(x) > |, we have 



oo oo 

\X 3 (x)\ < Ar l {2 U-\y)dy U\z){-- u{z))dz 



i v 
i y 



2 j r 2 (y)dy I <j) 2 {z){w{z)-\)dz 







< A^{\j + I}. (£.19) 



Combining (B.17)-(B.19), we derive 



A n 1 l 

J r < {2/ n+J _ 1 (0) K + (K - -)}( J + /(!)) + - J + 7. (£.20) 



From (A. 68) we have 

9 VV 3 4 



_ 24 , fWg 7. 
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Using (PT7p -(3.18) and ( gj§ ), we find for any m > 1, 

< r 1 ? = f ; < B - 21 > 

therefore 

An 1 1 

< R = {^ + K -lY J + I M) + - 4 J + L ( R22 ) 

Since for g sufficiently large, J, 1(1) and / are all 0( " g ) << 1, so is i? << 1. The Lemma 
is proved. 

From (|3.2| ) and (B.21), when n = 1, Aj satisfies 

A) = max - < / 1+i (0) + A(0) < 2K/£ 1 . (5.23) 
Thus, (B.22) implies 

A™ < RA^ 1 < R 2 A]- 2 < ■ ■ < R n ~ l A) < 2R n ~ 1 K/S 1 . (5.24) 

Consequently, given any e > 0, there exists an N such that (2R N ~ 1 K/£ 1 ) < e, and therefore 

Af = max \f N+j (x) - f N (x)\ < e. (5.25) 

for all j. Hence, by using Cauchy's test of convergence we establish the existence of lim f n (x) 
for all x. From ( |2.28|) , it follows that lim E n also exists. Theorem 3.6 is then proved. 
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Appendix C 



This appendix is divided into several sections. The proof of Theorem 5.2 is given in Section 
C.l. In order to set the bound for /C, defined in (|5.60| ), we need first the bounds for several 
other relevant functions and integrals. The bounds for ?p_, ip'_ and x' are given in Section 
C.2, the bound for X is in Section C.3, the bounds for J + and J~ in Section C.4 and the 
bound for Ai in Section C.5. 



C.l. Proof of Theorem 5.2 

As we shell see, the bounds a 7 and /3 7 in ( j5.19|) -( |5~2~0l) 



\2q 4 
7 < 4W— e"3f(l + a 



7> 4gJ*jL e -&(l-h 



can be written as 



where L = 0(-\iiy/g) is bounded by ( |4.45|) , and 



(C.l) 
(C.2) 

(C.3) 



1 - p y = [(1 + a M )(l + «tv) + ^J-e-^J+l- 1 



(C.4) 



with J + given by ( |4.27|) and (|4.29| ) and o^, a N are in turn bounds for 



M 



J (j)- + 2 {x)dx < ^ e ^(! + «m) 



(C.5) 



and 



which are given by 



CXJ , 

iV = j (pl(x)dx < \j^( 1 + a N) 



l + a M = (l- e~ 2 ^)e^(l + — f + 6g(e~*^ - e"^) 



(C.6) 



{C.l) 
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and 

l + a N = (l- -^)- 2 (l + A) + (1 + e -f). /i J_ e -** + * " 2 * (C.8) 
4^,3 48# V 7T g-Vb 2^2^ 

Therefore, for (7 sufficiently large 

a M = 0(g~ : 5), a N = 0(g~3) 

a 7 = 0(-ln^) and /3 7 = 0(</~s). (C.9) 

From these expressions, we anticipate that some of their derivations can be somewhat tedious. 
From CT, 

00 00 

7 - x = 2 / v; 2 (y)<fo / ^ 2 (z)^. (CIO) 



Take any I between and 1. We have 

1-1 1+1 

7 _1 >2 / Vv 2 (y)^ / ^ + (z)dz. (an) 



i-z 



By using ([0>p, (|0]) and (|42f3|), we see that 

^ m+[z) ~ my) (i + y?(o)(i + .) 2 • (ai2) 

Since 5 (y) = k 3 - y + | > -y + ~, 



,..7 -v :>, - y 1 3 

1-Z 1-Z 



J e 29S °( y) dy > e% 9 J e~ 29y dy = ^ 9 (1 - e" 2 ^ 1 " ). (C.13) 


Next, set z = 1 =F £ and write 5 (^) = £ 2 =F |£ 3 - We find 

Te^W . _ f 2(4 + acosh(|^ 3 )+8esinh(|^ 3 ) 

i z (i+-) 2 / (4 -a 2 e ^ 

> / se-^de > 7#(1 - "4= V 29 ' 2 )- (CM) 
J 2 4V2^ 1 ^hTg I } K J 

Combining (C.11)-(C.14) and using ( |4.44|) , we derive 

-1 . (1 - L Y ±g [X n _ -2j(l-I)ui _ 1 L-2 9 «2 



1 ^ ^ efg — (1 - e~ 2 ^)(l - -7= T^" 29 ' )• 
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which, for I = |, gives the upper bound (C.l) and (C.3) for 7. 

To derive the lower bound for 7, we observe that for y < z, in place of (C.12), 

00 

Substituting the inequality into (C.10) and decomposing the y-integration into two parts: J dy 

1 

1 

and J dy, we derive 



1 00 1 00 

7- 1 < J + + 2 J <P~\y)dy J <j>\{z)dz < J + + 2 J (j)+ 2 {y)dy J <p 2 + (z)dz (C.15) 
y 00 



where, according to ( |4.27|) and (|4.29| ), 

J+ = 2 / 0;%)^ | < — ln(e + 2e^2Vg). (CU6) 



1 2/ 



Using the definitions M and iV given by (C.5)-(C6), we see that the inequality (C.15) can be 
written as 

7- 1 < J + + 2MN or 7 > (J+ + 2MiV)~ 1 . (<7.17) 



Lemma 



and 



where 



(z) N<J^-(l + a N ), (C.18) 

(zz) N >^(l-p N ), {C.19) 
(m) M < ^et 9 (l + « M ) (C.20) 



H M> l e ^(l-/3 M ), (C.21) 



! + Qw = (1 - JjJ-d + A) + (1 + ^ + ^e-* (C.22) 
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l + a M = (l- e^ 9 )e^{\ + JLf + 6(e"3^ - e"^) (C.24) 

1 - /3 M = 1 - e ~ 29 (C.25) 
so that for 5 sufficiently large, a^v, /3jv, «m and /9m are all small, with 

a N = 0(g-*), p N = 0(-^-e- 29 ), (C.26) 

a M = 0(g-±) and (3 M = 0(e~ 29 ) (C.27) 

Proof of the lemma 

The lower bound for iV can be most easily derived by setting x = 1 =f£, as in (C.14), and 

00 

similarly observing that N = J <p\dx satisfies 



2 1 , — 

N > [ 7 ^e- 29So{x) dx > 1 2e~ 29e (% > J— (I - -== e' 29 ) (C.28) 

which gives (C.19) and (C.23). 

For the upper bound of N, we decompose N into a sum of three terms: 

N = N 1 + N 2 + N 3 , (C.29) 

where 



2 

1+/ 

N 2 = 

1-1 

and 

1-1 



N 1 = J <f) 2 + (x)dx, (C.30) 
2 

1+/ 

y (cm) 



1+i 

Note that, setting x = 1 + £, we obtain 



N 3 = ( J dx + J dx)<t> 2 + {x). (C.32) 



00 00 
iVi = / (^V^-S^df < / e-^df < ^e^. (C.33) 
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I 

For N 2 , the integral can be written as f(4> 2 + (x + ) + <j) 2 + {x_))d^ with x± = 1 =F £, as ' n 

o 



(C.14). Since 



4 2„ t 3 . 4 _2„ t 3 s o„f2 /■ 8 



by choosing I to satisfy 



we have for 



and therefore 



Likewise, from (C.32) 

) 2 
N 3 <8j e- 2 ^ 2 cosM-^K 



i 



1 | e- 2 ^ 2+ I^ 3 (l + e-t^ 3 )^ 
i 

< 4(1 + e-t 9 ' 3 ) | e-!^ 

oo 

< 4(1 + e-i^ 3 ) | e-t 9 « 2 ^ 



(C.34) 



\gi z = (^)\ (C-35) 



£ < ' = ^ (^-36) 
(1 " f r 2 < (1 - \ (C37) 
cosh(^ 3 ) < 1 + (^e 3 ) 2 , (C38) 



3 



^<2(i-^r 2 / e - 2 ^ 2 (i + ^ 2 m 

OO , 

< ^ - ^rj^a + = (i - ^)- 2 /|d + (a*) 



2gl g 1 ' 6 
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in which, because of (C.35), e~t s ' 2 = e' 3 ^ 3 . Combining (C.33), (C.39) and (C.40), we 
derive the upper bound for N, given by (C.18) and (C.22). 
Next, we write the integral for M, defined by (C.5), as 

a 1 

M={Jdx + J dx)(^^) 2 e- 29X+1 3 9xA+ ^. (C.41) 

a 

Choose 

a = 1/y/g (CA2) 

so that ga 3 = a = and 

a 1 

M < ( i + a )2e§3a 3 + |9 j e -2 9Xdx + y ( i + 5)2 e -| fl x+| flda . 

a 

< (l±^)2 e |^+| 3 l (1 _ e -2 9a) + /■ e -! 9 x+! flda . 

2 2 (? i/ 



a 



= — el 9 {(l - e" 2 ^) e 575(l + — f + 6(e"3^ - e"^)}, (C.43) 

which gives (C.20) and (C.24). 

For the lower bound, we have from (C.41) 

M > J - e' 29X+ l 9 dx = — et fl (l - e- 29 ), (C.44) 
o 9 

which gives (C.21) and (C.25) and completes the proof for the lemma. 

From (C.17)-(C.18) and (C.20), we derive the lower bound for 7, determined by 

7" 1 < ^\f^- e$ 9 [ (1 + «m)(1 + a N ) + 4«p e-t 9 J + ], (C.45) 
which leads to (C.4). 

C.2 Bounds for ^ and x' 
(i) Let 

1 

M(x) = J r+(y)dy (C.46) 



58 



and write 
with 

S(x) = S (x) + - In ±±f = i(x - l) 2 (x + 2) + - In ±±f (C.48) 
g z 6 g A 

We shall first establish: for 

0<x<l-#-\ (C.49) 
< (C7.50) 

Proof For < x < y < 1, define 

3(y) = . (C7.51) 

1 — x 

One can readily show that 

S(y)<?(y). 

Therefore, 

i 



Aiw< | e ^) ds< _^ (e ^,_ 1)< _^^ w , 



which gives (C.50). The restriction of x < 1 — g^ 1 in (C.49) is to avoid the uninteresting 
complication of S(x) having a zero at x = 1 — |<7 -1 - For 1 > x > 1 — since S'(x) < g~ 2 , 
we have </>+ 2 < e 2 ^ 9 , and therefore on account of (C.46) 

Mix) < - e 2 ' 9 = - (l + 0(~) ). (C.52) 
9 9 



(ii) Next, we turn to ip-(x). From ( |5.5|) , 

oo oo oo oo 

^ = 2 7 | ^;%)cfy | < 27/ 0;%)*/ / 

+ a; 3/ x y 

since ^ + = and < /+(?/) for z > y. Thus, for x > 1, 

t < 7 j +) (C.53) 
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where 7 S 4c/y^f e~3 9 and J+ < ± ln(e + 2e^huj) are bounded by (C.1)-(C2) and ( g2g) . 
For a; < 1 

< + 2 7 / 0;%)*/ < 7^+ + 2iV.M(x)), (<7.54) 

x y 



x 



where JV = / ^(^tfe is bounded by (C.18)-(C.19) and M{x) by (C.50) and (C.52). Hence, 



for 1 > x > 0, 

ij„(x)ij + (x) < ^ 2 + (x)(J + + 2NM(x)) = 0(e~3 9 ). ((7.55) 
It follows then, for all x > 0, 

^_(x)^+(x) < 0(e~&). ((7.56) 

At x = 0, on account of (|5.6|), 

V»-(0)V+(0) = ^.(0)=4/i(0)e-k 



Since ^ + increases from 0(e sf) at x = to 0(1) at x = 1, ^_ decreases from 0(e 3#) at 
x = to 0(e~t 9 ) at x = 1. 
Recall that, from (2.6)-(2.7), 

0_(x)0 + (x) = (_?_)V** (C.57) 
1 + x 

which, at x = 0, is 

0_(O)0+(O) = 4e"^. (C.58) 

Its derivative is 



1 +x 

at x = 0, the two terms on the left hand side are 



0_(O)0' + (O) = 4((7-l)e-^ 
and 

0+(O)#_(O) = -% + l)e-S'. (<7.59) 

Furthermore, for < x < 1 



< 0, 
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and therefore the derivative of the difference <p + — 0_ satisfies 

> <f>' = -gS'4> + . (C.60) 



U Ml ^ ±1 



The above expressions (C.55)-(C.58) describe the similarity between and ip_(x)ip + (x) 

As we shall see, the derivatives ip'+i 4>- ar| d x' = 4>'+ ~ i''- a ' so satisfy inequalities similar to 
(C.59)-(C60). 

(iii) In this subsection, we shall establish that, for < x < 1 — \ —, 

\ J V 9 

X \x)>i,' + {x)(l-0{e-y)). (CM) 
Proof From ( |4.43| ), ip + = f + 4>+ > <f> + . By using ( |575| ) we see that 

oo oo 

= ^-^'+ -27 / ip 2 + (z)dz < - 2 7 y" </>+(z)cfe, (<7.62) 

a; x 

which, at x = 0, becomes 

^+(0)^(0) < V- W;(0) - 2 7 iV, (C63) 

where, on account of (C.19) and (C.23), 

00 , — 
W = /^ W(fa >^-(l-_ e -»). (C.64) 



Because of ( ggg) and (C.47)-(C.48), 



V>+ /+ 0+ 0+ 1+3 

Multiplying (C.56), 

^+ < l^l(J+ + 2NM), (C.QQ) 
by ip'+/^+ and using (C.65), we find 

< 7</ + V>+'0+ + 27^^^+^+ < lJ+tp+tp+ - ZjNMiplgS'. 
Substituting this expression into (C.62), we derive 

M'_ < lJ+i/>+i//+ - 27(AT + gNMvplS'), (C.67) 
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where 

M(x) = I (j) 2 Az)dz. {cm) 



At x = 0, 

M(0)=M (cm) 

and 

jV(0) = N. (C.70) 

where M and A/" are given by (C.5)-(C6). 
On account of (C.65) and ([TP]), 

^+(0)^(0) < Vj(0)(p - 1) = /J(0)4(s - 1) e-t 9 < (1 - L)- 2 4(^ - 1) e"^ - 4<? e"K 

(an) 

where L = 0{-ln^g). From (C.1)-(C2), we see that 7 = 0(e), where e = e~ts as before. 
Therefore, 

7 J+^(0)=O(e 2 ) and 7 J+^ + (0)^(0) = 0(e 2 ) (C7.72) 
can both be neglected in (C.66) and (C.67). Because 



2 7 AA(0) = 2 7 iV > 2 • 4^/^(1 - ^)J—(1 - M e^ 9 = 8g e^ 9 (C.73) 



and 



-^(0)^(0)^(0) < ±- e*'(l + a M )(l ~ L)-H e^ 9 {g - 1) = \ , (C74) 
it follows then 

V>-(0)V>+(0) < 7 J+^+(0)^(0) - 2 1 NMgi) 2 + (0)S"(0) = 4# e~i 9 ((7.75) 

and 

V + (0)V_(0) < 7 J + V+(0)VV(0) - 2 T JV(1 + </MV£(0)S"(0)) 

= -8g e~l g {l -h = -4g e~l 9 (C.76) 

similar to (C.59). 

Next, consider the region 



< x < 1 - J^l (C.77) 
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so that the integral (C.68) for Af(x) already contains the Gaussian peak of (j)\(z) in its 
integrand. Hence, in the region (C.77), by following the same reasoning given in (C.28) and 
writing x = 1 — £, we find 



Af(x) = J <p 2 + (y)dy> J <p 2 + (y)dy>2je' 2 ^dy> ^(1 



1 1 



\-i \-i o ' a v 

For a; within the range (C.77), the smallest value of £ is y^^. Therefore 

In what follows, we shall first show 

C(x) = tp + (x)ip'_(x) — 7J + ?/> + (o;)-?//^(a;) 



- e 



)■ 



{c.n 



((7.79) 



is negative in the region (C.77). 

From (C.67) by using (C.78) for A/"(x), (C.50) for M(x) and ( p~9|) for ip+(x), we see that 



C{x) < -2 7 iV (1 - f3 x ){l + a^)- 1 + 



i , g( l-x)S'( x) 
2^(x) 



For g sufficiently large we can approximate 1 — /3jj , 1 + ax and f+(x) by 1. Hence, 

g(l-x)S'(xy 



C{x) < C(x) 



1 



2gS(x) 



in which 



g(l-x)S'(x) l (l-x) (l-x)(l-x 2 ) 3(1 + x) 

2^(x) 2 S^x) ol ' 2(l-x) 2 f(2 + x) 2(2 + x) 

Therefore for g sufficiently large, and within the region (C.77), 

C[x) = ip + (x)ip'_(x) — •yJ + ip + (x)ip' + (x) < 0. 



(c.so) 



(CM) 



(C.82) 



(C.83) 



We now turn to the derivative of 



X{x) = il>+(x) - ip-ix). 



Since 



Mx)x'(x) = MWAx) - M*W-{x) = ^+(x)#(x)(l - 7 J + ) - C(x). ((7.84) 
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On account of (C.83), 

tP+(x)x'(x) > MxWM£ ~ 1 J +) = M^+W ~ 0(e)), (C.85) 

which completes the proof of (C.61). As we shall see, this result greatly simplifies the deriva- 
tion of the bound for J. 



i y 

2/ 



C.3 Bound for X 
We begin with ( |5.58|) , 



J = 2 J x ~\y)dy J xWMz)**- 



As we shall establish, 

J<2(^) 1/4 (l + «x), (C86) 
where 1 + otj = 1 + 0(g~^) is given by 

1 + a x = 1(1 - L)- 2 + i[l - |(37r^)-4] • {l - l(37rp)-i - |(37r)^-f[l - l^)-*]- 1 

-^(S^l^fa-L)- 1 }" 1 (C7.87) 
and L = (^(Mn^/p) is bounded by (4.45). 

Proof Divide the integration range into two regions: 1. c < y < 1 and 2. < y < c. 
Correspondingly, 

J = l 1 +X 2 , (C.88) 

where 

1 2/ 
c 



and 



c 



l 2 = 2 I x -\y)dy I X (z)Mz)dz. (C.90) 

choose 

r: = 1 - , 
2 



c = 1 - 1(3^)^3 . (<7.91) 



64 



Hence, in region 1, 



MV) > Mc) = / + (c)rf e-WJW = (e-i(*) 1/2 ); 



1 + c 

whereas, according to (C.56) and for g sufficiently large 

and therefore can be neglected. Correspondingly 

i y 

Jx = 2 1 ^; 2 (y)dy / (C.93) 

c 

In region 2, on account of (C.85) 

X \z)>^' + (z)(l-0(e)). 
y 

Neglecting 0(e) and using x 2 (y) — 2 J x( z )x'( z )dz we have 

o 

y y 
c r Ix(z)^+(z)dz c J x (z)ip+(z)dz 

J 2 = 2 / dy^- y < 2 / dy\ < c/A, (CM) 

J o 2j X (z)x'(z)dz t 2$x(zW + (z)dz 



where 

tb' (z) 

A = min-±H f o r < z < c. (CM) 

To give an explicit form for both X x and X 2 , we assume first g very large so that only terms 
with the dominant power in g will be kept. Thus 

i y 
l L = 2j <b- 2 (y)dy J 4> 2 + (z)dz. 
o o 

From (A.62), 

<f>+ 2 (y)<f> 2 + (z) = (^) 2 e -MSo(y)S (z)) < 4 e -\ g (y-z? ^ (ag6) 
1 ~\~ z 

which leads to 

Ji<2^y(l-c). (C.97) 
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Likewise, in region 2, 

^+(z) 

Thus, 



-^)>^(l-c 2 ). 



X 2 < -r: ^ = t-t: r. (C.98) 

#(1 - c 2 ) 2#(1 - c) 

The optimal c is 1 — |(37rg)~3, which yields the leading form of J given by (C.86). The 
inclusion of non-leading order terms is tedious, as will be given below. 

For our purpose, we only need to establish a finite upper bound for X. For g large, it can 
be readily verified that the neglect of 0(e) terms is fully justified. Thus, from (C.93) and 
(C.97), we obtain 

Ji < 2J y (1 - c)/*(0) < 2W y (1 - c)(l - L)- 2 (C.99) 

where L = 0(^\n^/g) is bounded by (4.45). 

Likewise, (C.94) and (C.95) are valid after the neglect of 0(e) terms. For J 2 , we observe 
that 

From (4.55) and using f + (0)/f + (z) < (1 - L) _1 , we have 

TTT^-lS*- 1 - 1 )- 1 - (ai01) 

U{ z ) 2 f 9 

In addition, 

^M=g(l- Z 2)-J— (C.102) 

which at z = c = 1 — ^(3irg)~* gives 

^±l=g(tocg)-i [i--(3^)-4--(37r)i 5 "(l- 7 (3^)-i) 

= 2^(l-c)(l + 0(^)) . (C.103) 

Hence, 

J * < WT^f + a ^ (ai04) 
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where 

1 + a T2 = [1 - l -{^ g )-\][l - 1(3tt<7)-* - ^(37r)^-t[l - \{^g\ 

-^ViCl-L)" 1 ]" 1 . 
Combining (C.99) with (C.104)-(C.105), we have for 1 - c = §(37r#)-3 

„ . 3tt . i , . 



where 1 + a x is given by (C.87). 



C.4 Bounds for J+ and J" - 
From (|5.59| ) and ( |5.78|) , we have 



oo oo 

2/.AJ.. /\,2/ 



j + = 2 J x~ 2 (y)dy J r{z)dz 

i y 



and 



! 2/ 

2/ \ t . /\,2/ 



.7- = 2 / / • 



For x > 1, according to (C.53) 



^ < 7 J + = 0(e) . 



Thus, on account of ( |4.29| ) and neglecting 0(e), we have 



J + = J+< 7^1n(e + 2eJ2vr^ ). 
2g v 



For J7 , because x < V'+i we have 



J- < 1 < 2(4)^(1 + «x) • 



Combining together and J" , we derive 



J = J + + J~< — ln(e + 2 ev /2^) + 2(4)^(1 + «x) 
2# v # d 
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C.5 Bound for Ai 

Set n — 1 in ( |5.43| ). Since k = 1 and, according to ( |5.10| ), v\ = —7^+, we find 



/ Xtp+dx I X^+dx 



Ax =7^ 



< T- 

j 



((7.110) 



where, for convenience, we choose 



/ = 1 



'In g 



so that (i) for x > I, as can be readily seen from (C.53) and (C.56), 



ip + (x) 



0(e) 



(am) 



(C.112) 



and (ii) x = I is already outside the Gaussian peak of <p 2 + (x), which is near x = 1 with a 
width 0(-^=). Using (C.112) and neglecting 0(e), we can approximate the denominator in 
the righthand side of (C.110) as 



00 00 00 
X 2 dx = J ip+dx > J cf) 2 + dx 

1 1 1 



(C.113) 



For the numerator in (C.110), we note that 



00 00 

2j„. f2fn\ f j2 



therefore, 

00 

J X^+dx < J ip 2 + dx < fl(0) J (p 2 + dx 
00 
From (C.110) we derive, after neglecting the 0(e) correction, 



(C.114) 



/ <t> 2 +dx 



A 1<7 / 2 (0) 



I (p+dx 



((7.115) 



From (C.6) and (C.78), we have 



7T 



4>Ax)dx < J — (1 + ajv) 
2g 
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OO , 



I 



2g\ g 2 V2F[^< 
and therefore, for large g, A x = 7, with an upper bound given by 



<(l-L)- 2 (l + « JV )(l-^-=_ , (C.116) 
where L = 0^±j@) and a N = 0{g~*) are given by ( ggg ) and (C.8). 

Appendix D 

D.l In this section we examine the first iterated solution f\ and its associated energy S\, 
given by ( |2.31| ) and (2.32). From ( [2.31] ), we see that the derivative of f x is given by 



00 

f[(x) = 2^\x) J <f> 2 (y)(w(y) - £ 1 )dy, (D.l) 

X 

where because of (2.32), 

/i(0) = 0. (D.2) 

The energy E\ and the functions fi(x) and f[(x) are expressed in terms of integrals which are 
well defined for any g > 0. However, it is of interest to expand f[ and £% in power series of 

00 1 11 

/l = Z ~°Wl = _a 2 + — °3 H 

1 # m r 

and (D.3) 

00 1 11 

£1 = Z ~^ e m+i = ei + -e 2 + -^e 3 H . 

V 9 9 

As we shall see, these expansions give only asymptotic series. 

For such power series expansion, we can neglect terms 0(e~^ 9 ) and approximate (2.8) and 
(D.l) by 

<f>{x) = </>+(x) (DA) 



and 



00 

f[(x) = 20; 2 (x) / <j>\(y)(u(y) - SJdy, (D.5) 
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where (f> + and u are given by (2.6) and ( |2.11|) . Thus, the power series expansions of /i( 



x 



and Si are identical to those of f+,i(x) and S +i i, introduced in ( |4.15| ) and (|4.16| ) when n — 1, 
provided x is restricted to x > 0. We write 

fi(x) = U,i(x) and 

where f +: i(x) and satisfy 

which, on account of </> + = exp(— (y-So + Si)), leads to 

</£o/+,i = -f+ A - S[f' +)1 - S +tl + u. (D.6) 
At x = 0, according to ( |4.19| ), we have 

similar to (D.2). 

Since / +) i(x) and £ +j i have the same formal power series expansion as fi(x) and Si, we 
may write in accordance with (D.3) 

oo -y oo j 

i y o y 

Equating the terms proportional to y-~ n on both sides of (D.6) and recalling that 

uix) = (1 + x)~ 2 , we derive 

sX = (TT^- ei ' {D - 9) 

S'ov's = - - e2 (D. 10) 



So^m+1 = o^m - TTZ a 'm ~ e ™' ( D - U ) 



and for m > 2 

1 

2 Lm 1 + 
Since 

S^ = x 2 - 1 (D.12) 



is zero at x = 1, it follows from (D.9) that 



ei = - (£.13) 
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and 

Likewise, from (D.10) 
and 

a'=- 



^-h^Th +1) - (m4) 



e 2 = ^ (D.15) 



2 2 2 
5 ) 3 + 9 ) 2 + 9 ) + 9 



2 &y l + x' 

It is of interest to compare these results with the exact power series expansion of £ and f(x) 
derived in II. The above e±, e 2 , <j 2 and a' 3 are the same as the —Ei, —E 2 , —S' 2 and — S' 3 
of II, since, as can be readily verified, they satisfy the same equations. Note that at x = 0, 
f[(0) = f' +)1 {0) = 0, whereas a' 2 {0) = -§, a' 3 {0) = -^f , etc. Hence, the asymptotic 
expansion (D.3) for f[ or f' +1 fails completely at x = 0. 
In general, for m > 1, we may write 

i <-) 2m— 1 r> 



04r, 



Theorem D.l 



and 



2m- 1 

a L (m + l)= Yl (Z + 4)a,(m) (D.18) 

/ = max (0, L-2) 



2 4 

Proof From (D.17) and S[ = (1 + x)~\ we find 



1 „ „, , 1 2 ^ , w 2 



j<C-i " Vm + i = 24^+2 E «K^)( TT ^) i+3 (/ + 4). (D.20) 
Since its value at x = 1 equals e m+1 , it follows then 

i i 2m— 1 r> 

2Ci-^ + i-W = ^- E a,(m)[(— ) ,+s -l](Z + 4). (D.21) 

Using 

( 2 y+3 _ j = ( 2 _ 2 /+2 2 , +1 x] ^ D 22) 
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we obtain 

1 r\ 2m- 1 <") r\ 

and derive (D.18) and (D.19). 

It is of interest to construct the following pyramid structure of cfy(m): 

ai(l) «o(l) m — 1 

a 3 (2) a 2 (2) «i(2) a (2) m = 2 

a 5 (3) a 4 (3) a 3 (3) a 2 (3) «i(3) a (3) m = 3 (D.24) 

By using (D.18), we see that each row cti(m+l) can be obtained from the row cti(m) above. 
For example, using «i(l) = a (l) = 1, we have 

a 3 (2) = ai(l).(l + 4)=5, 
a 2 (2) = a 1 (l)-5 + a (l)-4 = 9, 

etc. The resulting elements in the pyramid assume the values 

11 m = 1 

5 9 9 9 m = 2 

35 89 134 170 170 170 m = 3 

315 1027 1965 2985 3835 4515 4515 4515 m = 4 (D.25) 

Correspondingly, (D.19) gives (D.13) and (D.15), 

170 85 4515 . 

e 3 = ^r = ^ 5 e4 = ^ir' etc - ( D - 26 ) 

We note that in this analysis e m is determined by setting the right hand side of (D.ll) to be 
zero at x — 1; otherwise, the corresponding cr' m+1 would diverge at the same x — 1. On the 
other hand, according to (D.l) and (D.3), the sum J2 a ' m +i( x )l9 m at a; = 1 is 

oo 

/;,!(!) = 2 / <^)[-£+,i + (7-r) 2 ]^ P-27) 
which is finite for any given set of e m which makes the sum 

OO J 

= E 7^ e ™+i (-D.28) 

o y 
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finite. 

D.2 In this section, we begin with the integral expression for f' +1 (x): 

oo 

f' +>1 (x) = 20; 2 (x) J <Pl(y)(u(y) - S + , x )dy, (D.29) 

X 

where 



OC 



/ <j)+(x)u(x)dx 

= ^ (£-30) 

JT <j) 2 + (x)dx 



in accordance with ( |4 . 15|) and (|4.16| ) for n — 1. Their expansions (D.8) lead to the equations 
(D.9)-(D.ll) for their expansion coefficients a' m and e m . Recall that in the above section, 
these equations are integrated from x — oo to x — 1. As remarked at the end of last section, 
we require each integrated solution a' m+1 (x) to be well-behaved at x = 1. This leads to 

1 9 85 .„ , 

ei = -, e 2 = ^, e3 = ^9' etc ' {D.31) 

with the corresponding cr' m+l (x) given by (D.14), (D.16) and (D.17). On the other hand, the 
function f' +1 (x), defined by the integral (D.29), is well-behaved at x = 1 for any finite 
as can be seen by the following explicit integral derived from (D.29): 

oo 

4,(1) = 2 / e-h(y-my + 2 ){ 2_ n _ s + 1 2 
7 1 + 2/ 1 + 2/ 



In the following, we shall show that unless the expansion £ +1 — ~^ e m+i has the same 

9 

coefficients e m+ i given by (D.31), the above f' +1 (l) does not have an asymptotic expansion 
in g^ 1 ; instead, it would have an expansion in . 
Introducing 

y-l = z, 

we first express (D.32) as an integral over z: 

oo 

f' +A (l) = 2jdz e- 29z2 [-a(z)S + , 1 +(3(z)] p.33) 



where 



a{z) = e-^ 3 ( -A_) 2 {DM) 
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and 



Next, we decompose a(z) and f3(z) as sums of even and odd functions of z: 

a(z) = cx+(z) + a-(z) 

and 



where 



and 



(3(z)=(3 + (z)+(34z), 
a±(z) = ~[a(z) ±a(—z)] 



(3 ± {z) = -[(3{z)±(3{-z)\. 



We then expand a + (z) and f3 + (z) as power series of z 2 . The results are 

a+(z) 



and 



+(i + 0(±))</V + (1 + 0(±))</V° + (| + 0(i)),V 2 + 



1 r, 5 9 ,4 35 1. 4 ,2 5 1 1 NN 9 fi 
+ ^ + ( 3 + 2^ +( 9 + 3^ + 0( ? )) ^ 

+( J? + O(i))^ 8 + + 0(i)),V° + (| + 0(i)),V 2 + 
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Correspondingly, the integrals 



and 



oo 

/i + = J dz e~ 2gz2 a + (z) 
o 

oo 

u+ = J dz e- 2gz2 f3 + (z) 



are given by 



and 



v , = 



1 /7T 



1 + 



35 1 ,105 7 2 -5 2 , 1 
+ 



^ 2\/ 2^^ 1 3-25 ^ ' v 2 9 ' 3 2 -2 n V 

1 HTr 89 1 395 59 • 9 • 7 • 5 • 3 1 1, 
8V2^^ + 3~25 ^ + ( 3~27 + 3 5 • 2 11 V + V' 



74 



We now introduce a new criterion, requiring the coefficient of each power g ^ n+ l) 



in 



oo 

^+ ■ Yl — e m+i + (DA6) 
o ST 

to be zero. This leads to the same set (D.31) derived before. The resulting power series 
expansion of (D.33) now contains only integer power of g^ 1 . 

Lastly, we observe that (including 0(e~% 9 )) is given by (D.30). To derive its asymp- 
totic expansion in g" 1 , we can neglect 0(e~l 9 ) and write 

<ft(x) = e~ 29{x - l)2 a {DAI) 



with 

1 + X 



oo oo 



In terms of x — 1 = z, (D.30) becomes 

oo 

S +A = J e- 29z2 p(z)dz j J e- 29Z ' 2 a(z)dz , (DA8) 
-l -l 

where a(z) and (3(z) are given by (D.34) and (D.35). In the decomposition (D.36) and (D.37), 
the odd functions a^(z) and f3-(z) contribute only 0(e~^ 9 ) to which do not appear in 



the asymptotic expansion of £+.,1. The yir/2g factor in (D.44) and (D.45) cancel each other 
in the quotient (D.48). The result is the same asymptotic expansion given by (D.23)-(D.26), 
as expected. 

In II we compute the coefficients E m in the asymptotic expansion of E ev , given by (^2|): 

OO 

E ev = g + J2 —E m+ i (£-49) 
o 9 

with 

1 9 89 

Ei = — A , E 2 = --, E 3 = - ¥ , etc. (£>.50) 

In Table 1, we list the ratios of \E m \/e m for several selected m from 1 to 100. 
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Table 1 

m \E m \/e m m \E m \/e. 



1 


1 


40 


1.5712 


2 


1 


50 


1.5692 


3 


1.0471 


60 


1.5668 


4 


1.1103 


70 


1.5643 


10 


1.4357 


80 


1.5620 


20 


1.5589 


90 


1.5599 


30 


1.5708 


100 


1.5579 



For m » 1, the E m is [10] 
and, from Table 1, 



E m ~--{\)™m\ {DM) 

71 O 



l(Vm! . (D.52) 

7T O 
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